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Abstract. A new, globally regular model describing a static, non spherical gravitating 
object in General Relativity is presented. The model is composed by a vacuum Weyl-Levi- 
Civita special field - the so called gamma metric - generated by a regular static distribution 
of mass-energy. Standard requirements of physical reasonableness such as, energy, match- 
ing and regularity conditions are satisfied. The model is used as a toy in investigating 
various issues related to the directional behavior of naked singularities in static spacetimes 
and the blackhole (Schwarschild) limit. 



1. Introduction 



Analytical modelling of astrophysical gravitating systems in General Relativity is a 
poorly developed field, especially due to the difficulties occurring in the field equations 
when the - obviously unrealistic - hypothesis of spherical symmetry is abandoned. As a 
' matter of fact we do not have any analytic model which could resemble a realistic (i.e. 

CO . rotating, axially symmetric) star This problem has an interesting theoretical counterpart 

' which is unsolved as well. It is the problem of finding sources for the rotating blackhole 

solution under conditions of physical reasonableness. Besides its obvious intrinsic inter- 
est, finding this kind of solutions would allow physicist to investigate those interesting 
(•~^ ', phenomena which can be expected to occur "at the boundary" between the singularity the- 

orems and the blackhole no-hair theorem. For instance, we do not know anything about 
those processes which should produce the progressive "cuts of the hairs" required to bring 
a non-vacuum, non singular rotating and collapsing configuration to settle down to a Kerr 
state before (naked) singularity formation. Of course the same question arises in the prob- 
lem of collapse of a non-rotating source of a Schwarzschild black hole, which according 
to the Israel theorem H] is the only static asymptotically flat vacuum solution to the Ein- 
stein's equations with a regular horizon. In order to delve deeper into these questions, we 
present here an interior solution which satisfies all physically relevant requirements and 
matches smoothly on the boundary surface to one of the static axially symmetric vacuum 
solutions 1 2 1 known as the gamma metric 1 3 1 . This metric belongs to the family of Weyl's 
solutions, and is continuously linked to the Schwarzschild space-time through one of its 
parameters. The motivation for this choice stems from the fact that the exterior gamma 
metric corresponds to a solution of the Laplace equation (in cylindrical coordinates) with 
the same singularity structure as the Schwarzschild solution (a line segment). In this sense 
the gamma metric appears as the "natural" generalization of Schwarzschild space-time 
to the axisymmetric case. This new interior adds to the few other known sources of the 
gamma metric |4| and would allow to study the behavior of very compact self-gravitating 
systems, which according to the Israel theorem are very sensitive to small fluctuations of 
spherical symmetry. 

2. Static, axially symmetric gravitational fields and the gamma metric 

The static, axially symmetric metric can be written, in full generality, as 
(2.1) ds^ = -x^dt^ + a'^lidx^f + (dx^)^] + P^dip^ 
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with Of = alx^jx"^),/] — (3{x^ ,x'^),x — xi^^i^'^)- The components of the Ricci tensor 
satisfy to 

(2.2) i?* + i?^ - 

where A/ = g"^iyi + g"^iyi is the flat Laplace operator in the fictious two-dimensional 
space parameterized by cartesian coordinates x^ and x"^. 

Thus, if we are considering vacuum gravitational fields (so that Ricci is zero) or if the 
matter source satisfies to + — = then one can make use of the so called Weyl 
gauge. This consists in changing variables (a:^, x^) {g, z) where g := (3x is harmonic 
and z is its harmonic conjugate. It follows that the metric can be written as follows: 

(2.3) ds^ = -e^^df + e2''-2A [(^^)2 ^ (^^^2] ^ g2^-2X^^2 

where 

(2.4) a = e''-\x = e\ 

In the present section we shall consider only metrics satisfying the Weyl gauge (the physi- 
cal content of this condition will be discussed in next section). The field equations can, in 
this case, be written as follows: 

(2.5) kT^ = 2gX,gg + 2gX^z^ + 2X,g ~ gX^^ - gX^^ - gv.gg - g^/^z^ 

(2.6) KT^g = g{Xl-Xl)-v^g 

(2.7) nT', = -Q{>^^^a^>^%) + ^^e 

(2.8) - 2gX,gX^,-u^, 

(2.9) kT:^ = -g{X^^g + X%+u,gg + y,,,) 
where we put 

K-.^^^ge""-^^. 

It is worth recalling the structure of the vacuum solutions. For vanishing energy-momentum 
tensor, we get 

(2.10) AA = 0, 

(2.11) = e(Afe - X%),v^, = 2eA,eA,,, 



(2.12) Az/ + A^^ + A^^ = 0, 

where A is the flat Laplace operator on the fictitious two-dimensional space parameterized 
by the cylindrical coordinates g and z (i.e. A/ — f^gg + f^^z + Once a solution of 
the Laplace equation has been chosen, the second equation becomes an exact differential 
and allows the calculation of 1/ by a quadrature: 

(2.13) / g[{X\~ Xl)dg + 2X^gX,,dz]+c 



r 



with r any open curve in the {g, z) plane, and it can be verified that the third equation 
identically holds. There is, therefore, a one to one correspondence between solutions of 
the Laplace equation on flat two dimensional space and vacuum, static, axially symmetric 
spacetimes. Asymptotic flatness requires vanishing of A at space infinity, and the general 
solution of Laplace equation is thus of the form |5 1: 

00 

(2.14) A = ^-^P„(cos^) 

n=0 
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where 

(2.15) R = Vg^ + z'\ cosi/j = ^ 

and Pn{cosip) are Legendre's polynomials. It then follows: 
n^f.\ „ ^ (n + l)(fc + l) a„afc 

n,k=0 

The real constants a„ are the so-called Weyl moments. The Weyl moments are not the mul- 
tipole moments measured by an observer at space infinity (6| . For instance, the Schwarzsch- 
ild solution of mass m has the following Weyl moments: 

(2.17) a2„ = -- — — ,a2ri+i=0 

2n + 1 

and thus it is not the monopole-Weyl solution. It is however easy to relate the a„ (Weyl) 
and Mn (multipole) moments for instance it can be shown that 

(2.18) Mo = M = -ao, 



(2.19) M2 = g = 02 - ^4. 

In what follows we shall concentrate on that particular metric corresponding to the solu- 
tion of the Laplace equation for a source of constant density 7 uniformly distributed on a 
segment of length 2m on the z axis centered at z = L3J: 

(2.20) A = l,jR. + R--^rn 



2 \R+ + R^+2m 

(2.21) . . 2^^JiR,+R--2m)(R,+R^+2rn) 

^ ^ 2 \ AR+R^ 



where i?+ = + {z + rn)^, R_ = ^J + [z — m)^. This metric is at best visualized 
in the so-called Erez-Rosen coordinates {p, z) {r, {}): 

(2.22) = (r^ - 2mr) sin^ z (r - m) cos 1} 
in which: 

7 / 2m\ 

(2.23) A = ^ In 1 - 



2 V r J 
(2.24) V = —In' 



2 yi^^ + ^sin'dj 
so that 

(2.25) ds^ ^ A'^'-'^-^E^-'^'dr^ + A'^'-^S^-'^'dt^^ + sin^ dA^-'^dip^ ~ A'^dt^ 
where 

(2.26) A= l ,S=1 + — sin^ ■& 

\ r J \ J. J.Z 

This metric is usually called gamma metric . It has many interesting properties which we 
are going to discuss below. 
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3. Physical properties of the gamma metric 

The gamma metric coincides with the Schwarzschild metric of mass m when 7 = 1. 
One of the advantages of the Erez-Rosen frame is that in this limit the gamma reduces to 
Schwarzschild in Schwarzschild coordinates. The mass and the quadrupole moment of the 
gamma metric are given by Isj: 

(3.1) M = 7m, 



(3.2) 



In the 7 = limit the metric becomes flat, whilst the Curzon solution (i.e. the Weyl mono- 
pole solution) is obtained in the limit 7 ^ 00, m ^ with 7m = const \9\. 
As we have seen the gamma metric is originated - formally - by the solution of the Laplace 
equation for a line bar of constant density. However this a rather formal identification, 
and we are interested here in properties of the gamma metric as a vacuum field observed at 
space infinity. First of all, one can study the surface of revolution analogue to Schwarzschild 
spheres BlOI . 

To an r = rt surface in Schwarzschild pertains an area 



(3.3) 



A, 



Jo 



2tt 



r1 sin 'dd'ddip — Aiirl 



Calculating the same area for r = rf, surfaces in the gamma case, we can classify such sur- 
faces as oblate or prolate with respect to the Schwarzshild sphere. Obviously, this property 
of the revolution surfaces is also the expected "shape" of static axially symmetric sources 
of the metric. 

The 2-dimensional metric on the slices r = ri,,t = const is 



(3.4) 



dn' = r?AZ 



2 • 2 

r,, sm 



where a subscript denotes evaluation at r = r?,. One has, to lowest order in e = 7 — 1, 



(3.5) 



detg^ 4 sin^^Al' E^"'^ = sini?(l - elnSb). 



From which follows: 



(3.6) 



I - -e 
2 



siniJlnShdt? 



Since In Ef, < we get > As for e > (7 > 1); in such cases the slices are oblate 
(respectively prolate for 7 < 1). If the metric has to be interpreted has the vacuum field 
generated by a symmetric object then the physically interesting case is obviously that of 
oblate sources (e.g. galaxies). Stability has, for instance, to be expected only for such 
objects. This result is in agreement with the physical shape since the expansion of the 
newtonian potential is given by l;6J: 



500 



(3.7) 



-1 



= -1 



^ + ^P.(cos.) + a(l 



7TO 
r 



1 7771"^ 



(Y - 1)P2(C0S79) +0 — 



so that the newtonian potential of the source is that of a disk if 7 > 1. 
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4. The nature of the singularity of the gamma metric 

In the present section we investigate on the nature of the singularity of the gamma 
metric, extending and completing previous results by Virbhadra 1.1 1.1 . The Kretschmann 
scalar K = Rabcd.R"'''"^ is 

(4.1) KM- 16r.VA^(''''^) 



rl0^272-27+25]3-272 

where 

N{r,'d) := sin^ ??[3r7i7(7^ + l)(m - r) + 7^(4to^ - 6mr + 3r^) + 
(4.2) +m^(7'' + 1)] + 3r{mj + m - rf{r - 2m). 

For 7 = 1 we get the Schwarzschild result K — ^^"^ , which is obviously regular at 
r = 2m. In what follows we consider only the case 7^1. The r = 2m surface is then 
a true singularity of the spacetime and the no-hair theorem assures that it will be a naked 
singularity in the sense that it is not covered by an event horizon. 

The diverging behavior of X in ?■ = 2m obviously depends on For = one has 
E = A and therefore: 

_ 48mV(m7 + m-r)2 
^^■^) ^8A4-27 ' 

this function diverges at r = 2m only if 7 < 2. 

If ,9 = ^0 ^ X(r,i9o) ~ A-'^'+T-i. Since 7^ - 7 + 1 > V7 7^ 0, the Kretschmann 
scalar diverges in r = 2m if the singularity is approached on any plane different from the 
polar one, for any value of 7. It is therefore interesting to investigate the visibility of such 
singularity in dependence of the direction in which it is approached. 
Consider a null "radial" (i) = const) geodesic. A faraway observer can receive a signal 
from the r — 2m surface if 



+0 



(4.4) lim / dt = lim / A^ ^ ^S^dr < 



2m+e ^ J2m+e 



00. 



On any plane different from the polar one (-d ^ 0) there exists a positive constant C such 
that 

/■^ 2 2 /""^ 2 

(4.5) lim/ a( "' ~V~^ ')Yr^dr <limC A^ ^ "'"'"' ^dr < 00 

since the integral on the right converges for any 7 7^ 1 (if 7 = 1 divergence of the integral 
obviously signals the presence of the event horizon). 
If however ?? = 0: 

pR ,2_ _ _ 2 fR 

(4.6) lim/ A(' 2"' ')s-^dr = lim / A-^dr 

and the integral converges if and only if 7 < 1. The situation is summarized in the table 1. 
We conclude that the singularity is naked independently on the direction of observation if 
7 < 1 while it exhibits directional nakedness for 7 > 1, being not visible, in this case, on 
the polar plane. 





7 < 1 


7 > 1 


t9 = 


visible 


not visible 


^^0 


visible 


visible 



Table 1 . Directional behavior of the singularity of the gamma metric 
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5. Sources of the gamma metric 

To construct physically viable sources of the gamma in order to investigate their be- 
havior "at the boundary of the no-hair theorem" one needs to chose a matter model and 
then integrate the coupled Einstein-nmatter field equations. However, this proved by now 
to be impossible in the non-spherical case (also in the case of the Schwarzschild field few 
exact interior solutions are known). We thus need a simplifying assumption which how- 
ever should restrain as little as possible the physical significance of the results. There is, in 
fact, the famous argument known as "Synge trick" which proposes to "guess" the metric 
from general considerations and then to evaluate the matter content from Einstein equa- 
tions; however, usually such "Synge solutions" have nothing to do with physics. The idea 
we follow here traces back to work by Hernandez (TT\ and it is based on deformations of 
physically valid spherically symmetric solutions. We interpret here these deformations as 
actual changes in the shape of the objects, i.e. we use e as the deformation parameter. 
As far as we are aware the unique known examples of interior solution for the gamma are 
that given by Stewart et al. 1 4 1 . One of this solution is based on the Schwarzschild constant 
density interior solution as a "seed", but the resulting density is not strictly decreasing from 
the center at higher orders and therefore the solution cannot be stable with respect to radial 
perturbations. The second one is based on the Adler interior solution. 
In order to construct a new interior solution we take here a different point of view. Rather 
than searching for a solution generated by a "seed" to lowest order, we consider a "varia- 
tion of the mass" of the gamma metric, namely, that family of interior solutions which is 
obtained by (a priori) any inhomogeneous distribution of mass which generates the vacuum 
gamma field at a fixed boundary. It is obvious, that not all the mass distributions will be 
physically viable, and indeed we shall see that the space of allowed functions is severely 
restricted by physical conditions. 

It is well known that an interior solution generating the Schwarzschild field can be obtained 
from the vacuum metric putting m ~ IJ-{r). The Einstein tensor then becomes: 



(5.1) 






(5.2) 


G\ 




(5.3) 


Gl 


At) 


(5.4) 


Gl 


f/{r) 
r 



So that the matching, regularity and energy conditions 



(5.5) 






m 


(5.6) 


lim ''^'J 







(5.7) 


4 


> 





(5.8) 


j,2 J, 


> 






can easily be satisfied, for instance, by a matter distribution of the form 
a9, ,(r)^i.iJ..-=(l-|^). 

Note that the energy condition Gq — G\ > holds identically for all r. 

We now apply the same reasoning to the gamma. Considering 7 = 1 + e to first order in e 
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the vacuum metric becomes: 



(5.10) 
(5.11) 
(5.12) 
(5.13) 



goo 

911 

922 
933 



-A - eAln A 

A-i + eA-i(lnA-21nE) 
+er^ (In A - 21nS) 
sin^ 1? — er^ sin^ In A 



So that taking the Schwarzschild interior solution as the order zero interior solution is 
possible to build an interior of the form 



(5.14) 



(5.15) 



500 = - 1 - 



Mr) 



911 = 



= 1- 



r 

Mr) 



-e 1- 



1 



Mr) 
r 

Mr) 



F{r) 



{F{r)+G{r,^)) 



(5.16) 322 = r'' + er''iF{r)+Gir,§)) 

(5.17) 333 = sin^ 1? - er^ sin^ ^?ii^(r) 

Matching conditions require 
(5.18) 



(5.19) 
(5.20) 

(5.21) 
(5.22) 

(5.23) 



^J'irb) 
f^'irt) 

F{rb) 
F'irt,) 
G{n,i9) 



m 


In ( 1 



2m 



'b 
-2 In 

4m 



2m \ 

rb J 
2m 
rb 
2m 

n 



m 



sin^ § ] = A{d) 



while regularity conditions for r 

(5.24) 

(5.25) 



require 

MO 

G{r,d) 



+ 








Now the energy conditions have to be imposed. It turns out that a suitable solution which 
satisfies the weak energy conditions can be found only for e > 0, i.e. for oblate sources. 
The first energy condition, namely Gq > 0, is satisfied automatically for all r ^ if the 
corresponding condition for Schwarzschild is satisfied. In r = the order zero condition 
vanishes and so [GgJ^^rb ^ must be imposed. From this follows: 



(5.26) G"{rb,'&) < 1- 



2m 

n 



4m 
"3" 



>G^^{n,'&) 



1 



m 

rb 



G'in,^) 



In the same way conditions Gq — G| > and Gg — G3 > are satisfied automatically for 
all r ^ and so only [Gg - Gi]r=o > and [Gg - G|]r=o > must be imposed. The 
first is however satisfied by means of the regularity condition for G(r, 1?) while the second 
leads to 



(5.27) 



lim 



2E1^+F"{r) ) = c>0. 
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Finally the condition that must be satisfied for every r and -d is 

G",-Gl = (^l-M!:))(^V(r)-iG"(r,,?)) + 

(5.28) (GMr,^) + —^G,{r,i))j > 0. 

It is easy to check that the latter is satisfied in r = and r = r?, if the previous conditions 
are satisfied and so, by continuity it must be satisfied in a neighborhood of those two points. 
Taking an appropriate function G{r, ■&) so that G"{r, iS) and G^^{r, ??) + ^^G^{r, are 

small compared to for all i} satisfies the condition everywhere. Explicit examples of 
choices of F{r) and G{r, ■&) satisfying the conditions are given in the appendix. However 
we show in the figures the behavior of the energy density and the critical energy condition 
for one of such examples describing a compact object with m = IAMq and r?, = lOKm. 




(a) (b) 

Figure 1. In figure are shown (a) the condition for the energy density 
Tg > and (b) the condition Tg — > (in the adimensional unit 
obtained multiplying by r^) for a fixed value of e = in dependence 
o/'i? and the adimensional radial variable 



6. Discussion and conclusions 

General Relativity is in principle the "final theory" for the description of macroscopic 
gravitating objects. It is therefore odd enough that after nearly a century from its invention 
we do not have even one sound model of a star within the theory. In fact no solution of 
stationary axially symmetric Einstein field equation in matter is known, which is able to 
represent a finite non singular rotating object. Even in the simpler case of static, but not 
spherical, gravitating objects very little is known. Of course this case is the non rotating 
limit of the previous one, since it is only in the case of the perfect fluid that the static object 
must be spherical, and rotating stars are not to be expected to be made out of perfect fluids. 
For instance it is very well known that anisotropy plays a relevant role 1131 in the structure 
of compact objects such as white dwarfs and neutron stars. 

So motivated we investigated here a simple static but not spherically symmetric model 
of a compact object, showing the existence of physically valid solutions of Einstein field 
equations describing it. 

This model should allow for studying the dynamic behavior of the compact object close 
to the surface ri, — 2m. This may be done, either by considering the dynamic equation 
just after the system leaves the equilibrium 1 14| or by calculating the active gravitational 
(Tolman) mass as the system leaves the equilibrium L15 1. In either case the bifurcation 
between the exactly spherically symmetric case and the gamma solution, should appear 
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even for very small values of 7, if ^ is sufficiently close to i. 

However, observe that specific constraints appear in our model on the value of ^. Indeed, 
the central radial pressure is given by 



P(0) - Eo+e -Eo hi 1 



2 m 



F'in) 



n 



Eorl 



C 



1 

8^ 



which is consistent with the small departure from spherical symmetry approximation if 



--Bnlnl 1 



F'in) (I 



n 



1 



C (\ 



-E^ri - — + 3 - — 



1 



6 



Stt 



< E^ 



Since both ln(l - and = 2™ (^1 _ diverge as ^ i, the pressure 

blows up in that limit. We must therefore consider an upper limit in the surface gravita- 
tional potential for which the model is acceptable. The same reasoning can be applied to 
the "arbitrary" positive constant C. In fact once the scale given by e is fixed the value of C 
must not exceed a limit value for which the approximation holds. 
Given the fact that ^ — ^irEor-^ maybe it is better to say that once we fix the value of 
(and so Eq) the possible values of e become fixed as well. So that for — - 



we must 



have e — > 0, and the more we depart from spherical symmetry the less the value allowed 
for — becomes. 

Finally we would like to mention that work aimed to extend these results to slow rotation 
scenario is in progress. 



Appendix A. Example 
The Einstein tensor for the interior metric l5"l 415 . 1 7l results : 



(A.l) 
(A.2) 



(A.3) 



(A.4) 



= 2 



1 

M'(r) 



fi{r) 



-G) 



r 



G' 



1 

2^ 



-2^{F + G) 



2r 



G'- 



2r 



G{){) 1 
"2^ " 2 

^ Mr) 
r 

cos 'd G{, 
sin-i? 2r2 

^ Mr) 



1 



Mr) 



G"- 



-2;(^'W- 



{t + Gj : — + 

Mr) 



r 

Kr) 



1(1 



r 
fi{r) 



{2F' - G') 

{2F' + G') 

{2F' + G') 
Mr) 



2r2 



(^4F' + G')--{ 1 



r 

Mr) 



{2F" + G" 



{2F' + G') 



A simple choice of F{r) satisfying all conditions is 



(A.5) 



1 F'in) 2 

Fir) = r 

2 n 



1 



2 r 



^G^[l~-- 

6 n 



D 



with G an arbitrary positive constant and D given by the condition for F(r) at the bound- 
ary. 

Furthermore considering G'(r, d) of the form 



(A.6) 



G(r,?9) = 5(r)A(z9) + /i(r)B(i?) 
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with g{r) > so that g{rb) = 1 and g'{ri,) ~ and h{r) < so that h{ri,) ~ and 
h'{rb) = 1 satisfies the matching conditions. The regularity requires: 

(A.7) lim 4) = lim = 0. 

r^o r 1 — »o r 

The condition — G\ > is then satisfied if g{r), g"{r), h{r) and h"{r) are small 
enough with respect to F(r). Which is always possible due to the degree of freedom given 
by the choice of the positive constant C. The condition [GQ]r=rb > finally impUes 



n J \rl rl \ n J n 

which is easily satisfied. 

If we consider an object with m ~ 1.4i\/Q and rt — IQKm (i.e. a neutron star) in ge- 
ometrized units we get — = ^ and Eq — (which corresponds to 2.6 • 10^^-^) then 

we have a wide variety of choices for the functions F{r) and G{r, d) and the central pres- 
sure doesn't diverge for every small value of e so that the reasonable choice e = leads 
to a central pressure ^^j^ times greater than the corresponding pressure for the spherically 
symmetric case. 
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